Quantum-entanglement aspects of polaron systems 
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We describe quantum entanglement inherent to the polaron ground states of coupled electron- 
phonon (or, more generally, particle-phonon) systems based on a model comprising both local 
(Holstein-type) and nonlocal (Peierls-type) coupling. We study this model using a variational 
method supplemented by the exact numerical diagonalization on a system of finite size. By way 
of subsequent numerical diagonalization of the reduced density matrix, we determine the particle- 
phonon entanglement as given by the von Neumann and linear entropies. Our results are strongly in- 
dicative of the intimate relationship between the particle localization/delocalization and the particle- 
phonon entanglement. In particular, we find a compelling evidence for the existence of a nonanalyt- 
icity in the entanglement entropies with respect to the Peierls-coupling strength. The occurrence of 
such nonanalyticity - not accompanied by an actual quantum phase transition - reinforces analogous 
conclusion drawn in several recent studies of entanglement in the realm of quantum-dissipative sys- 
tems. In addition, we demonstrate that the entanglement entropies saturate inside the self-trapped 
region where the small-polaron states are nearly maximally mixed. 

PACS numbers: 71.38.Ht, 03.67.Mn 



I. INTRODUCTION 

In recent years entangled quantum systemsii^ have gar- 
nered interest as a resource for quantum information 
processing."^ In addition, a great deal of research effort 
has been expended towards clarifying the role of entan- 
glement in the (zero-temperature) quantum phase tran- 
sitions of many-particle systemsi^>^i^ Perhaps the most 
compelling, however, is the need to elucidate its possible 
bearing on the macroscopic properties of physical sys- 
temsiii^ Namely, while the entanglement entropies are 
simply related to the many-body density matrix of the 
system, they bear no a priori relation to any observ- 
able physical quantity. Attempts to associate a physical 
meaning with entanglement were made, for example, in 
mesoscopic physics: schemes have been proposed for de- 
tectioii^ and even measurement of entanglement by ex- 
ploiting its emergent relation to the quantum noise i^^d^ 

Central to all the above developments is the problem 
of quantifying entanglement in diverse physical systems, 
such as quantum spin chainS j^^d^'^*^ interacting bosons 
and/or fermionS j^^d^d'^d^ quantum-dissipativ e^^'^° and 
disordered systemSf^i to name but a fewi^^ In this re- 
gard, one of the areas of condensed-matter physics whose 
quantum-entanglement aspects have heretofore received 
only scanty consideration is that of the polaron problem. 
Ever since its inception by Landau and Pekarj^S, the po- 
laron concep t^^'^^1^^ a quantum particle interacting with 
a bosonic environment - has played an immensely im- 
portant role in theoretical studies of coupled electron- 
or exciton-phonon (henceforth e-ph) systems. What 
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is more, this truly ubiquitous concept is lately finding 
resurgence in seemingly unrelated physical situations, the 
realm of ultracold atoms being a case in point 

The main body of polaron-related work is focussed on 
the study of a single electron interacting with the har- 
monic lattice vibrations through a short-range, non-polar 
potential that is linear in the lattice displacements and 
describes the dependence of the electronic on-site ener- 
gies on the lattice degrees of freedom. The traditional 
starting point in describing such interaction, dubbed lo- 
cal e-ph coupling, is the paradigmatic molecular-crystal 
model due to Holstein.— As is being amply appreciated 
latelyr^i^i^ however, local coupling is not the only type 
of short-range e-ph interaction relevant in realistic sys- 
tems: nonlocal (off-diagonal) e-ph coupling accounts for 
the phonon-modulation of the electronic hopping inte- 
grals and bears relevance to several classes of molecule- 
based systems. The most common form of nonlocal 
coupling is Peierls-type coupling j^^i'^^i^^ widely studied 
within the framework of the semi-classical Su-Schrieffer- 
Heeger (SSH) mode l^^'^^i"^^ that describes the anomalous 
transport properties of nonlinear excitations (solitons, 
polarons) along the quasi-one-dimensional polyacetilene 
chain. In an implicit way - becoming manifest by carry- 
ing out the Jordan- Wigner transformation - this coupling 
forms the basis of the XY spin-Peierls model^^i^ hence 
the name. Likewise, coupling to the breathing mode 
in cuprate superconductors is of Peierls' type and - as 
transpires from recent investigations based on general- 
izations of the t — J model ~ plays an important role in 
these systems, especially in the regime of weak dopingi^ 
Somewhat different forms of nonlocal coupling have been 
shown to be of relevance for charge transport in organic 
molecular crystals )^"'^^ carbon nanotubes^^ and DNA 
wires4^ It is worth of mention that a form of nonlocal 
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coupling has also been incorporated in a generic electron- 
boson coupling model recently proposed by Alvermann et 

Given the abiding interest in the polaron problem 
in condensed-matter physics, the quantum-entanglement 
aspects of this problem have so far not been given due 
attention. While the changeover from a small to a large 
polaron is known to have the nature of a smooth crossover 
(with no broken symmetry), rather than a phase transi- 
tion,*'^ it is still tempting to quantify it using the entan- 
glement measures. Entanglement in the one-dimensional 
Holstein model was studied by Zhao tt al^ The authors 
emphasized the relation between the self-trapping pro- 
cess and the quantum (hetero-) entanglement between 
the phonon subsystem and electronic excitation. Effects 
of Peierls-type coupling, however, are as yet totally unex- 
plored; given the wealth of intriguing implications of such 
interaction it appears interesting to elucidate the role it 
plays in particle localization as seen through entangle- 
ment measures. Additional motivation comes from re- 
cent investigations of entanglement in the realm of quan- 
tum dissipative systems One of the most important 
conclusions of these studies is that entanglement mea- 
sures can have non-analyticities away from any phase 
transition and that these non-analyticities are intimately 
related to the loss of coherence. To be more specific, 
Stauber and Guinesfc^ found a nonanalyticity at the tran- 
sition from underdamped to overdamped oscillations in 
the Ohmic case**^ of the spin-boson model. Moreover, 
this nonanalyticity proved to be even more pronounced 
than the one occurring at the actual localization phase 
transition. 

As a matter of fact, an evidence that non- analytic be- 
havior of entanglement-related quantities does not nec- 
essarily coincide with the quantum phase transitions had 
already been found before.* Besides, a study of local- 
izable entanglement in a gapped quantum spin systemic 
has showed that entanglement length diverges despite the 
fact that the correlation length remains finite, the latter 
indicating absence of a quantum phase transition. On 
the other hand, in disordered systems, for example, no 
such cases have been reported. It is therefore of interest 
to investigate whether the polaron problem - somewhat 
related to the spin-boson model, but with no phase tran- 
sitions taking place - may also defy the tenet whereby 
an occurrence of a nonanalyticity in the entanglement 
entropy is a telltale signature of a quantum phase tran- 
sition. 

In the present work, we study entanglement in po- 
laron systems. In order to address the problem from 
as general a viewpoint as possible, we start from a po- 
laron model that includes both local (Holstein-type) and 
nonlocal (Pcierls-typc) short-range c-ph coupling. Given 
that this Hamiltonian does not admit an exact solution in 
any physically-relevant limit, we analyze it using a varia- 
tional method supplemented by an exact diagonalization 
on a finite-size system. Proceeding in this way, we find 
that entanglement exhibits non-analyticities as a func- 



tion of Peierls-type coupling, which are not accompanied 
by a phase transition. We argue that the occurrence 
of such "accidental" nonanalyticities in the problem at 
hand is related to the loss of coherence, much like in the 
quantum-dissipative systems. 

The outline of the remainder of this paper is as fol- 
lows. In Sec. In] we present the model and notation to 
be used throughout. Sec. IIIII contains details of our 
variational approach: we first introduce our variational 
Ansatz, then provide details of analytical derivations 
needed to implement it, and finally describe our com- 
putational method of variational minimization. The fol- 
lowing Sec. IIVI contains essential details of the exact di- 
agonalization method. In Sec.|V|we introduce the entan- 
glement measures and lay out the method for calculating 
them. The obtained results are presented in Sec. IVIi ac- 
companied by a discussion of their salient features. We 
conclude, with some general remarks, in Sec. IVIII 

II. MODEL 

The system under study consists of an excess particle 
(electron, hole, exciton) interacting with harmonic lattice 
vibrations - dispersionless (Einstein) phonons - through 
a short-range interaction. As our starting point, we adopt 
a one-dimensional e-ph Hamiltonian obtained by dove- 
tailing the Peierls-type coupling term on the conventional 
Holstein Hamiltonian. (Restriction to a one-dimensional 
system is not a severe limitation given the short-range 
nature of the e-ph interactions discussed here.) The com- 
pact form of this extended Holstein model reads 

H = ^ e,({u}) a|a, -h ^ i^+i ,,({u}) {a\^^a, + h.c.) 

i i 

+ cjJ2blk, (1) 

i 

where a|(ai) creates (destroys) a particle at i-th site (at 
position z = 0, 1, . . . , — 1), &|(6i) creates (destroys) 
a dispersionless phonon with frequency lu at the same 
site. The effective on-site energy 

ei{{u}) = e + anUi , (2) 

and hopping integral 

ti+i^i{{u}) = -t + ap{ui+i - iii) , (3) 

depend on the lattice displacements Ui = {2Mcj)~^^^{bi + 
bl), where an and ap are the local (Holstein-type) and 
nonlocal (Peierls-type) coupling constants, respectively, 
and M is the mass of molecules in the underlying crystal. 
The bare on-site energy and hopping integral are denoted 
by e and t, respectively. For simplicity we take e = in 
what follows. More explicitly, the Hamiltonian can be 
written as 

H = He + ^ph + Hg + H^, (4) 
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where 

He = -t'^{al^-^^a^ + h.c.) , 

i 

i 



= (j)uj^{al^j^a^ + h.c. )(&!+! + k+i ~ H - 



(5) 
(6) 
(7) 

(8) 



with g = aH/V2Muj^ and cj) = ap/V2Muj^ being the 
dimensionless local and nonlocal coupling constants, re- 
spectively. The two important limiting cases of our model 
arc the Holstcin model {(p = 0) and the quantized version 
of the SSH model (5 = 0). 

The eigenstates of Hamiltonian ^ ought to be the 
good-quasi- momentum states, i.e., eigenstates of the to- 
tal crystal momentum operator 



(9) 



since the latter commutes with H. In the following, 
the eigenvalues of K will be labelled with k. By mak- 
ing use of the Born- von Karman periodic boundary con- 
ditions, the quasi-momenta in the first Brillouin zone 
are given by k = {2-k / a){ra/N) {a - the lattice spac- 
ing; m = 0, 1, . . . , iV — 1). For convenience, we express 
quasi-momenta in units of a^^, so that nRn — una —> kti. 



III. VARIATIONAL METHOD 
A. Choice of variational Ansatz 

While not being exactly-soluble, Hamiltonian in 
Eq. can be treated variationally; in the studies of cou- 
pled e-ph systems methods of this type have been shown 
to yield quantitatively trustworthy results that compare 
well with those obtained by the exact (numerical) diag- 
onalizations4^i^ An important class of such methods is 
furnished by Toyozawa's Ansatz stated and generaliza- 
tions thereof42ii^iS While these Ansatz states have been 
widely used in the studies of the Holstein model, they are 
also capable of describing systems with simultaneous lo- 
cal and nonlocal coupling. 

Omitting the most conventional form of Toyozawa's 
Ansatz state, we purposefully cast it in a way that ren- 
ders manifest its entangled nature: 
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|O)e0n^'(^'-J|O)ph, (10) 



where |0)c (|0)ph) is the electron (phonon) vacuum, and 
Di{a) = exp(a&[ — a*bi) is Glauber's displacement op- 
erator that creates the phonon coherent state |^;1„)ph = 



£'i(C;l„)|0)ph at site /. This (overcomplete) set of phonon 
coherent states captures the multi-phononic nature of the 
polaron ground state. Importantly, it is the dependence 
of the variational parameters 7m_„ on both m and n that 
renders Ansatz state entangled: the sums over these 
indices cannot be decoupled, implying that this state can- 
not be expressed as a separable (direct-product) state in 
the Hilbert space Ti. = He Tiph- 

Generally speaking, the use of variational methods in- 
variably involves the trade-off between flexibility of the 
variational wave- function (which increases with the grow- 
ing number of parameters) and numerical difficulty of 
finding reliably the global minimum of the ground-state- 
energy expectation value (complexity grows rapidly with 
parameter number; see Sec. IIII 51 for additional details). 
In this regard, the major drawback of Toyozawa's Ansatz 
is that it involves a large number - 27V for each k - of 
variational parameters. Alternative methods have been 
proposed that provide accurate results, while involving 
smaller number of parameters. Adopting this point of 
view, we seek the polaron eigenstates of Hamiltonian ^ 
in the form of translationally-invariant Bloch states 



\Mn)) 



(11) 



with "form- factors" |7/;„(n)) given by"* 



iK,7n 



|o). 
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))|0>ph- (12) 

The skew-Hermitian operators C/k;j (j — 0,±1,±2) are 
defined as 



J7«;,(n) = -^^(/.;, (q) e^«" 6,-h 



(13) 



(14) 



with 

^'"'^^'^^ " 1 + 2 {t/uj) /3«y[cos(K) - cos(k + q)] 

Using trial wave- functions pT|) . the lowest polaron 
band Eq<^{hi) can be obtained by minimizing the en- 
ergy expectation value over variational parameters = 
{aK;o, aK;±i, aK;±2; /3k;o,/3k;±i, /3k;±2; ^^(-S), . . . , ^^(S)}: 



Ecsi^) = min 



(V'«|g|V'.) 

(V'kIV'k) 



(15) 



[Note that for each k there are in total twenty variational 
parameters; this number is fixed, rather than being pro- 
portional to the system size as in Toyozawa's Ansatz. 
Importantly, among eleven parameters ^kU) only ten 
are independent because of the normalization condition 
on the trial wave- function in Eq. pT|) .] In the following, 
we will be particularly interested in the polaron ground- 
state energy Ei^s(k = 0) = Eq. The corresponding 
variationally-optimized ground-state wave-function will 
hereafter be denoted with | GS) . 
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B. Matrix elements and computational scheme 

Here we present the derivation of the expression for 
the energy expectation value, foUowed by the details of 
our numerical method for variational minimization. To 
facilitate further derivations we first note that 



exp 



( E 



i=-2 
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N 



where Dq{aq) = exp(a<j6j — a*bq) is Glauber's displace- 
ment operator that creates a coherent state of phonons 
with quasi- momentum q, and w^iq) is defined as 



It is straightforward to show that 



m—m' ' 



(17) 



(18) 



where 



mm' q 

X (e-'"'"'i«.(<?)-e^^'"<(g)) . (25) 
The above formulae are easily derived using identity 

gAgS ^ ^A+B^hM , (26) 

which holds if operators A and B satisfy condition 
[A, [A, B]] = [B, [A, 13]] = 0, as weU as identity 



(27) 



valid for an arbitrary analytic function f{b\ b) of bosonic 
operators. 

An important measure of the multi-phononic nature 
of the polaron ground state is the average number of 
phonons 



ph 



GSI J2 blh 



GS) , 



(28) 



-tqm 



N 



-iqm 



(19) 

By making use of the well-known expression^^ for the 
overlap of coherent states (a|/3) (a,/3 G C) 



= {0\b\a)b{P) |0) 



- p"*/3 +1/31 ) 



we obtain 



exp 



N ■ 



Other relevant matrix elements are given by 



(20) 



(^kI i?ph i^k) - ^ E *:k)i>.(™) 



(22) 



^e''(™-"') |i«.(g)p, (23) 



AT 

X ^ (e-'«"' u;, (q) + e^'^u): (g) ) , (24) 

9 



which can be obtained using Eq. 

Based on the expressions for | B |V'k) and (t/'kIi/'k;), 
we perform variational minimization in order to find the 
polaron ground state (k = 0) for a system with A^ = 32 
sites. Due to the large number (n — 20) of variational 
parameters involved, the energy expectation value is a 
function of these parameters with multiple local minima. 
Finding the global minimum thus constitutes a rather 
nontrivial numerical optimization problem. We perform 
this complex task using the multi-start-based global ran- 
dom search methods^ we first generate a large sample 
10^) of random points in the space of variational pa- 
rameters; we then select a smaller number (~ 20) of them 
that have the smallest values of the function to be mini- 
mized and perform local searches for minima \p(rc') com- 
putation] around each of these points: the one with the 
smallest energy is then adopted as the sought-after global 
minimum. The fidelity of this approach is corroborated 
by the stability of the final result for the global minimum 
upon varying the initial number of random points. 

Regarding the choice of our numerical method a re- 
mark is in order here. In the local-coupling-only case 
((/> — 0, i.e., the Holstein model), the system is ex- 
actly soluble in the "unphysical" limit of zero hopping 
{t — > 0) by the well-known Lang-Firsov canonical trans- 
formationi^l The natural way to proceed in finding the 
optimal variational parameters for finite t is then to start 
from the exact solution for t — Q and gradually change 
using optimal values of variational parameters obtained 
for given value of t as the initial guess for the next, 
slightly higher value. However, for finite the model 
is not exactly-soluble in any relevant limit. Therefore, 
the procedure just described does not carry over to the 
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4> case and one needs a careful sampling of the entire 
space of variational parameters, afforded by the multi- 
start-type methods, to reliably find the global minimum 
of the ground-state-energy expectation value. 



C. Application scope of our variational method 
and comparison with other methods 



IV. EXACT DIAGONALIZATION 

To supplement our variational method, we also per- 
form an exact numerical diagonalization of Hamiltonian 
(III) on a system of finite size. We study system of = 6 
sites, varying the maximal number of phonons between 
M = 8 and M — 10. The states in the truncated Hilbert 
space are given by 



While Ansatz state in Eqs. pT|) - p2)) was originally in- 
troduced for use in the Peierls-coupling-only case (quan- 
tized SSH model), we here demonstrate the that it can 
also be utilized for the local-coupling-only case (Holstein 
model), and accordingly for the case with simultaneous 
local and nonlocal couplings. The nearly-perfect agree- 
ment of our results in the ((> = case with another known 
variational approach, suggested by Cataudella et al.^^ , is 
illustrated in Fig. [T^,b. The latter was shown to agree 
well with the Global-Local method^^ implying the gen- 
eral agreement of all three methods. Therefore the trial 
states used here are generally applicable to short-ranged 
e-ph interactions with Einstein phonons. The agreement 
between the results obtained by the variational meth- 
ods of the present type and other approaches (density- 
matrix renormalization group methodpS quantum Monte 
Carlo^) is also well established.— 




FIG. 1: Comparison of the ground-state energies obtained 
using the variational method of Cataudella et aL^^ (circles) 
and the method of the present work (solid curve) in the (j) = 
case: t/u — 1.0 (a) and t/u = 2.0 (b). 



(29) 



are 



where n = (rio, rijv-i) and m = (TOo,...,m, 
the sitewise electron and phonon occupation numbers 
(Si "-2 — — ^^)- Coefficients Cn,m contain 

the information about the phonon content of state j^"). 
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FIG. 2: Phonon distribution sitewise in the polaron ground 
state for t/uj = 1.0, (j) = 1.3 (iV = 6, A/ = 10). The states 
m with the same total number of phonons m = Tn-i are 
indicated by the bins at the bottom portion of the plot. 



One of the crucial prerequisites for a successful applica- 
tion of the exact-diagonalization approach in the present 
context is a proper truncation of the (otherwise infinite- 
dimensional) phonon Hilbert space. In other words, the 
maximal total number of phonons on a lattice has to be 
large enough as to be capable to account for the phonon 
distribution in the polaron ground state in the strong- 
coupling regime. Typical phonon content (sitewise) of 
the polaron ground state obtained by exact diagonaliza- 
tion is illustrated in Fig. [21 where the n-th group of peaks 
represents the phonon distribution for the case when the 
electron is located at site n {n = 0, ...,N — 1). Within 
each group of peaks, coefficients Cn,m are given in the 
ascending order of the phonon occupation numbers. For 
weak coupling, phonon distribution pm ~ X^m l^n,mP 
(where m — pcaks at m = phonons, while 

in the strong-coupling regime it peaks at m = 3 or 4 
phonons. This a posteriori corroborates that our choice 
of the maximal number of phonons (between 8 and 10) 
was pertinent. 
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V. ENTANGLEMENT 
A. Entanglement entropies: generalities 

In order to make the present work self-contained, be- 
fore embarking on the calculation of entanglement in our 
coupled e-ph system we review the general prescription 
for characterization of bipartite quantum entanglement. 
We consider a composite quantum system that can be 
divided up into two parts A and B, where A denotes 
the subsystem of interest and B the environment whose 
details are unimportant. The Hilbert space of the full 
system has the form of a tensor product: Ti. — Ha Hb- 
In a pure state |^), the density matrix of the full system 
is given by 



MM 



(30) 



We then construct the reduced (marginal) density matrix 
Pa by tracing over the environmental degrees of freedom: 
Pa = Tr^ p. The von Neumann (entanglement) entropy, 
defined by 



S = -TrA(/OA In Pa 



(31) 



contains information about the quantum correlations 
present in the pure quantum state under study. It rep- 
resents the most widely used measure of bipartite quan- 
tum entanglement. Note that S = — TrA(/5A InpA) = 
— TrB(pB Inpg), where the reduced density matrix p^ is 
obtained by tracing over the degrees of freedom in subsys- 
tem A. The upper bound S'max — ln(£') of S, where D is 
dimensionality of the reduced density matrix, is reached 
when the reduced density matrix is maximally mixed. 
The so-called linear entropy is defined as 



S, = 1-Tva{pI) 



(32) 



and has the advantage (compared to the von Neumann 
entropy) of being easier to calculate. It is worth of men- 
tion that both the von Neumann and linear entropies are 
closely related to the quantum Renyi entropies^ ^gip) = 
(1 - (j)-iln(Trp«) {q > 0): S{p) is the q ^ 1 limit of 
Sq{p), while Sl{p) is simply related to S2{p)- The upper 
bound of Sl, reached for the maximally-mixed reduced 
density matrix, is given by ^^^max = 1 — D~^. 



B. Reduced density matrix from variational 
approach 

In accordance with general relation ([50]) . the density 
matrix corresponding to the state IV'k) on the tensor- 
product Hilbert space H = He ^ Hph is given by 



Pe-phin) 



(V'kIV'k) 



(33) 



The reduced particle (electron) density matrix is given 
by the partial trace over the phonon Hilbert space Hph- 



(34) 



/5c (k) = Trph[/Oc-ph(K)] 



Straightforward derivation, with the aid of Eq. (|19|) . 
yields 



nn' m,?n' 



X Z. 



m—m'-\-n' — n 



jK.{n-n') 



(35) 



|n)ee(n'| 



where \n)^ = |0)o is the state with electron at site n. 
The last equation, when combined with Eq. (jlSp . readily 
leads to the expression for the general matrix element of 
the reduced density matrix pc{n — Q) = Pc in the polaron 
ground state: 



1 E,„„>^Lo("^O^K^o(m) z^^,^^,^^ 



(36) 



The corresponding von Neumann entropy 

S = -Tro(poln/5c) (37) 
can be expressed as 

5 = -E^«l^^^' (38) 

i 

where {Xi | i = 0, 1, . . . , iV — 1} are the eigenvalues (A^ > 
, Ai = 1) of po. Based on our variational approach, 
the linear entropy 



(39) 



can readily be obtained in an analytical form using Eq. 
P6p : however, we here omit the ensuing cumbersome 
expression. In the exact-diagonalization approach, the 
reduced density matrix pc is obtained from the corre- 
sponding eigenvectors. The results illustrating depen- 
dence of the entanglement entropies on the e-ph coupling 
strengths are presented in the following section. 



VI. RESULTS AND DISCUSSION 

Our numerical results, obtained using variational ap- 
proach, correctly reproduce a continuous dependence of 
the polaron ground-state energy, displayed in Fig. [3^ for 
t/uj = 1.0, on local and nonlocal coupling strengths, re- 
flecting the well-known absence of phase transitions in 
coupled e-ph systemSi**- 

The entanglement entropies, depicted in Fig. ^p,c for 
t/uj — 1.0, both behave in a similar way: they in- 
crease with increasing coupling strengths, saturating and 
remaining essentially unchanged in the self-trapped re- 
gion where the particle becomes localized. The satu- 
ration values of the two entropies are close to those of 
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FIG. 3: Variational ground-state energy (a), the von Neu- 
mann entropy (b), and linear entropy (c) as functions of di- 
mensionless local (g) and nonlocal (</>) coupling strengths, for 
t/uj = 1.0 and iV = 32. 



the maximally-mixed density matrix: S'max — In(A^), 
•Si,, max = 1 " N''^. For example, for N = 32, S'max = 
3.46, S^,max = 0.97 (cf. Fig. ED. 

Importantly, for t/u; > 0.85 we find a strong nonana- 
lyticity in the dependence of entanglement entropies on 
the nonlocal coupling strength. This nonanalyticity has 
the character of a jump-discontinuity and is more pro- 
nounced in the case of von Neumann entropy (Fig. [5]d) 
than for the linear entropy (Fig. [3}:). It becomes more 
and more pronounced with increasing value of t/uj, i.e., 
upon approaching the adiabatic regime t ^ uj. In or- 
der to emphasize that discontinuous behavior sets in for 
sufficiently large value of the ratio t/ui, in Fig. 3^ we de- 
pict the von Neumann entropy for t/uj = 0.25 where the 
nonanalyticity does not exist at all, and for t/uj = 2.0 
(Fig. where the nonanalyticity is noticeably more 
pronounced than for t/uj = 1.0. 

To emphasize the appearance of a nonanalyticity in 
the entanglement entropies as a function of nonlocal cou- 
pling strength, we study the nonlocal-coupling-only case 
making comparisons between the results of the two ap- 




FIG. 4: The von Neumann entropy S, obtained variationally, 
as a function of dimensionless local (g) and nonlocal {(f)) cou- 
pling strengths. Displayed are results for t/uj — 0.25 (a) and 
t/uj = 2.0 (b), with N = 32. 



proaches (variational and exact diagonalization). Typical 
results are depicted in Fig. [5l Figure [5^ illustrates very 
good agreement between the two approaches as far as 
the ground-state energy is concerned. In Figs. [5)3, c the 
entanglement entropies are shown as obtained from two 
different variational calculations (for iV = 6 and N — 32) 
and two different exact-diagonalizations (with N = 6 and 
maximum M = 8 or 10 phonons used). The only sizeable 
discrepancy is in the behavior of 5, which clearly stems 
from the finite-size effects: while values of 5" correspond- 
ing to the variational calculation with iV = 32 deviate 
considerably from those of exact diagonalizations with 

= 6, the difference between the two approaches when 
applied to a system with the same number of sites is not 
very drastic. However, the most important feature of the 
obtained results, manifest in all the cases considered, is 
the non-analytic behavior of S and with respect to 0. 

Detailed analysis shows that the observed nonanalyt- 
icities occur for values of 4> at which the lowest en- 
ergy states of H{(j)) undergo avoided crossings. The 
crossings are avoided (rather than real ones) because 
[H, Hcf,] 7^ leading to a smooth dependence of the 
polaron ground-state energy on the coupling strength (j) 
(and, accordingly, the absence of a phase transition in 
the conventional sense of the term). There is, however, 
no general principle that would rule out the occurrence 
of nonanalyticities in the entanglement entropies at these 
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FIG. 5: Ground-state energy (a), the von Neumann entropy 
(b), and linear entropy (c) for tju) — 1.0 and g = 0, as cal- 
culated using the variational (var.) and exact diagonalization 
(e.d.) approaches. N is the number of sites, while M stands 
for the maximal number of phonons used for exact diagonal- 
ization. The arrows indicate the entropies for the maximally- 
mixed reduced density matrix. 



avoided-crossing points. 

The fact that nonanalyticities of the type discussed 
here do not exist in the case with local e-ph coupling - 
nor in the "statically" -disordered systems described by 
the Anderson model - and that for Peierls-type coupling 
they show up only when t is larger or of the same order as 
uj point to the possible importance of "dynamical disor- 
der" (retardation, i.e., nonlocality in time) effects, which 
are here bringing about the nonlocal particle-phonon cor- 
relations. Namely, when the relevant electron energy 
scale (set by the hopping integral t) becomes comparable 
to or larger than the characteristic phonon energy (w), 
the lattice deformation does not follow instantaneously 



the electron motion. Consequently, the phonon modes 
that are excited by the passage of the electron take a 
long time to relax. Therefore, a lattice deformation can 
be observed far away from the current position of the 
electron. As a consequence, the effects of retardation in 
the e-ph interaction become prominent. Such effects are 
known to be much more pronounced for the Peierls-type 
interaction than for the purely local Holstein-type inter- 
action, even when the effects of phonon dispersion are 
accounted for in the latter i22 This can be traced back to 
the fact that unlike local coupling, which is momentum- 
independent, the Peierls-type coupling depends strongly 
on both the electron and phonon momenta. More pre- 
cisely, in momentum-space this coupling reads 



+ b,) , (40) 



where 7(fc, q) is the e-ph interaction vertex function 

7(fc,g) = 1-;^= (sin(fc)-sin(/c + g)) . (41) 

In particular, at small phonon momenta 7 behaves as 

7(fc,<7)(xg (<Z--0), (42) 

which is a very strong momentum dependence, and dif- 
ferent than that of the Frohlich e-ph interaction.^^ Lo- 
calization of a "Peierls polaron" is therefore expected to 
have a much more dramatic impact on the nature of the 
accompanying e-ph correlations than that of a "Holstein 
polaron" . Additionally, the more pronounced character 
of the nonanalyticity observed for larger t/ui also appears 
to be in consistency with this argument; it demonstrates 
the increasing "inertia" of the more and more spatially- 
extended phonon cloud to electron's localization. 

In the light of our findings and those of Stauber and 
Guinea,— it is tempting to draw some parallels between 
the two models involved, or more specifically, between 
our model and 'Ohmic systems'. However, our polaron 
model - at least in its full form - does not seem to bear a 
direct relation to any of the known quantum-dissipative 
models, because of the intrinsically off-diagonal nature 
in the electron bilinear operators and the dispersionless 
character of phonons that it involves. [Besides, based 
on Eq. (UJ) we can infer that the Peierls '-coupling term 
with acoustic (rather than optical) phonons would be the 
most similar to the super-Ohmic systems with spectral 
density J (to) oc w^, even though we are here concerned 
with a one-dimensional system.] The link to these models 
appears much easier to establish for the local-coupling- 
only Holstein model (0 = 0): the two-site version of this 
model, discussed long time ago by Shore and Sander^^ 
represents a simplistic (single-mode) form of the spin- 
boson model. However, the nonanalyticities of entangle- 
ment entropies that we find in the present work occur 
in the polaron crossover regime, and are accompanied by 
the growth in the average number of phonons in the po- 
laron ground state [cf. Eq. (UHl)]- This can indeed be 



9 



considered as a physical situation analogous to the loss 
of coherence in the spin-boson-type models. 

A few remarks regarding our variational method are 
in order. While examples are known of artifacts^" in 
the variational approaches to coupled e-ph or quantum- 
dissipative systems - a prominent one being the failure 
of variational methods to predict a continuous transition 
in the sub-Ohmic case of the spin-boson model - in the 
case at hand such an approach correctly reproduces the 
smooth dependence of the polaron ground-state energy 
on both local and nonlocal coupling strengths (the po- 
laron crossover). The obtained nonanalyticities in the 
entanglement entropies, corroborated through the exact 
diagonalizations, represent a robust feature that under- 
scores the connection between entanglement and local- 
ization. 



VII. CONCLUSIONS 

In summary, we have investigated the quantum- 
entanglement aspects of polaron systems. As our point 
of departure, we have adopted a very general polaron 
model that includes both local (Holstein-type) and non- 
local (Peierls-type) particle-phonon coupling. We have 
studied this Hamiltonian using a sophisticated varia- 
tional approach supplemented by the exact (numerical) 
diagonalization on a finite-size system. We have estab- 
lished a close connection between the entanglement and 
phonon-induced localization. Our results make it trans- 
parent that the entanglement entropies constitute much 
more sensitive indicators of the change of polaron states 



than the ground-state energy. While the intuition as to 
the relationship between entanglement and localization 
has already been manifest from studies of other classes 
of systems^^ - most prominently the disordered ones - 
our findings add to it some other elements. 

As a salient feature, we have demonstrated that - 
above some threshold value for the ratio of the hop- 
ping integral and the phonon frequency - entanglement 
entropies exhibit a nonanalyticity as a function of the 
nonlocal (Peierls) coupling strength. This nonanalytic- 
ity is physically related to the loss of coherence and is 
not accompanied by a phase transition. In this sense, the 
present work reinforces the conclusions drawn in some re- 
cent studies of related quantum-dissipative systems, such 
as the spin-boson model. Furthermore, our findings un- 
derscore the fact that the nature of the phonon-induced 
localization in the presence of nonlocal particle-phonon 
interactions is different than that of the purely local in- 
teractions. This may have bearing not only on the solid- 
state systems exhibiting polaronic behavior, but possi- 
bly also on certain classes of cold-atom systems - with 
(in principle) tunable couplings - where phonons can be 
introduced in a controlled wayi^ The need to investi- 
gate the interplay between entanglement and the phonon- 
induced localization in other relevant models' is clearly 
compelling. 
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